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Abstract 



It is shown that practically any physically relevant random distribution of 
offset charges destroys the Kosterlitz-Thouless-Berezinskii charge-unbinding 
phase transition in two dimensional normal and superconducting granular ar- 
rays and films. The array conductance obeys the Arrenius dependence on 
temperature. Offset charge disorder decreases the effective Coulomb gap of 
the system and may account for recent experimental findings in two dimen- 
sional arrays of tunnel junctions. 

PACS numbers: 73.40 Gk, 73.40 Rw 
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Charging effects in tunnel junction arrays remain under intensive investigation during 
last several years. These effects become important at sufficiently low temperatures T < 
Eq = e 2 /2C where C is the characteristic capacitance of a tunnel junction between metallic 
islands and/or of such an island itself. Under this condition tunneling of a single electron 
from one island to another is essentially suppressed due to the Coulomb interaction and the 
system conductance G shows an activation behavior G(T) oc exp(— aEc/T), where a ~ 1. 
At T = electron tunneling is blocked completely and the system becomes insulating. 

The above arguments are quite general and apply to the most of mesoscopic granular 
arrays. It was pointed out by Mooij et al. that in two dimensional (2D) arrays a new 
interesting collective effect may take place. Provided the island capacitance to the ground 
Co is small as compared to the junction capacitance C electrostatic interaction of charges 
H] in such arrays logarithmically depends on the distance A between these charges up to 



A ~ JC/Cq, A is measured in the lattice spacing units. Thus in the limit Cq — > the 
system represents an example of a 2D Coulomb gas which exhibits a Kosterlitz-Thouless- 
Berezinskii (KTB) phase transition 0|| at a temperature Tktb of order Ec- Physically it 
implies that at T > Tktb there is a nonzero concentration of free charges in the system 
and its conductance remains finite G > 0, whereas at T < Tktb & U charges are bound in 
charge- ant icharge pairs and the system linear conductance drops to zero G — 0. Thus 2D 
arrays may become insulating not only at T = but also at finite T < T KTB ~ E c /4:TT |],@|. 

The above conclusion applies to both normal and superconducting granular arrays. In 
the latter case an elementary charge is that of a Cooper pair 2e (instead of e for normal 
arrays) and therefore Tktb for superconducting arrays is four times larger than that for 
normal arrays 

Several experiments were performed to study charging effects in both normal and 
superconducting 2D granular arrays and none of them indicated the presence of a KTB phase 
transition for charges. E.g. no specific KTB dependence G(T) ~ exp(— A/^/Tktb — T) in 
the vicinity of Tktb ~ Ec has been found. In contrast for a wide temperature region the ar- 
ray conductance was reported to follow a purely activation behavior G(T) oc exp(—aEc/T), 
where for most of the samples the parameter a varied between 0.23 and 0.27 []I],|6|-[8]] . These 
results might look somewhat surprising. Of course in real systems finite size as well as the 
selfcapacitance effects turn a charge-KTB transition into a crossover. Nevertheless for the 
sample parameters [l],P-|8| one estimates A ~ 10 4- 20 and thus such a crossover could be 
expected to be sufficiently sharp to be distinguished from a purely activation behavior. 

An important assumption made in [[| was that the island charge is quantized in units of e 
(or 2e for superconducting arrays). In real experimental situation, however, this assumption 
is frequently violated because of the random charges trapped in the substrate. These charges 
can polarize the islands and induce noninteger offset charges on them. Some consequences 
of this effect were previously studied numerically an d experimentally [|TT|,|I1|. 

In this Letter we argue that randomly distributed noninteger offset charges have a dra- 
matic impact on the behavior of 2D granular arrays and that the experimental results 
|]J6H8| acquire a natural physical explanation provided the effect of offset charges is taken 
into account. We demonstrate that already a very weak disorder completely washes out the 
charge-KTB transition in 2D arrays leading to an activation-type behavior of the system 
conductance G(T). We also show that in the presence of the offset charges the Coulomb 
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gap for electron tunneling across a half of junctions in the array decreases for any given 
direction of the applied current. At sufficiently low temperatures these junctions yield a 
dominating contribution to the system conductance and the effective Coulomb gap becomes 
smaller than that in the absence of the offset charges. 

Let us consider a square tunnel junction array with the junction capacitance C and 
assume that the selfcapacitance of the islands is very small Co — > 0. The partition function 
of this array can be written in terms of a path integral (see e.g. 
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where the collective variable 4> x {r) is linked to the voltage V x on the island x by means of a 
standard relation = 2eV x (r), Q x is the island charge, (3 = 1/T and the term S t [4>} describes 
electron tunneling between islands. In a superconducting array the variable </> x (t) represents 
the phase of the superconducting order parameter of the island whereas in normal arrays it 
is not more than a formal variable. The sum in (|5|) is taken over all integer winding numbers 
n x on each island x. 

In what follows we shall assume that quenched offset charges Q x are randomly distributed 
over the system, so that Q x = and 
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Depending on the physical situation different types of disorder can be considered. E.g. 
the values Q x on each island can be fixed by strong local potentials in the substrate being 
completely independent of each other. Then we have g(x — x') = g u 5(x — x') or, equivalently, 



9k — 9u, 



(3) 



gk is the Fourier component of g(x — x'). Alternatively one can assume that Coulomb 
interaction between offset charges on different islands dominates over the local potentials. 
In this case we have 



g k = 4 sin 2 (fci/2) + sin 2 (A; 2 /2) g c = 4A(k)g, 



(4) 



where for the model adopted here [15] g c = 2ir 2 T/Ec and the components ki 2 of the wave 
vector k are normalized by the inversed lattice spacing constant 1/a. We believe that 
practically any physically realistic situation can be described by a proper combination of @ 
and (|). 

Let us first assume that tunneling between islands is small and disregard the term St in 
the expression for the partition function (|l|). Then following |2[ we make a shift x (r) = 
<f) x o + AiTn x T/ (3 + 6 x (t) and perform Gaussian integration over the variables 8 x (t) with the 
boundary conditions 6^(0) = 9 x (/3) = 0. Then we get 
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where we denote K = 2tt 2 T/Ec- For the sake of definiteness let us first consider the case of 
a correlated disorder (|) . In order to average the free energy F over disordered offset charge 
configurations we make use of the replica trick and write 
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F = -T —Z m 

dm 
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Performing the standard transformation to the field theory [14] with the aid of the Poisson 
resummation formula we arrive at the expression 



Z m = J Dip exp | - 
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where (p is the m-component field. We introduced here the chemical potential hiy which 
unrenormalized value is set equal to zero. 

In order to proceed with the renormalization group (RG) analysis we come to a continuum 
limit (p x — (p x > —>■ aW(p(x) and introduce the minimum length scale £ ~ a. The idea of 
calculation is essentially the same as that of ||. The generalization concerns only the effect 
of disorder. Expanding the exponent in the expression for Z m in powers of y and g and 
successively integrating out the Fourier components f^, with q — 5q < q' < q one comes to 
larger and larger scales e l and arrives at the RG equations 
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Here we have to put g = g c . 

An analogous calculation can be carried out for the uncorrelated distribution of the 
offset charges (|^). In this case the last term in the square brackets of (|6|) has the form 
— (<7u/2) J2x{J2 Q n x) 2 - Repeating the same procedure we again arrive at the RG equations 
(§) with g = g u . 

It is interesting to point out that the scaling equations (§) have the same form as those 
derived earlier in Ref. [110] for a 2d XY model with random phase shifts. However an 



important physical difference is that in Ref. |16] the case of random dipols (or, equivalently, 
anticor related charges) was considered whereas here we are dealing with random offset 
charges Q x . Apparently at large distances random charges should be more efficient in 
screening the charge- ant icharge logarithmic interaction than random dipoles. And indeed 
- although the disorder term in our equations (^) has the same form as in Ref. |TB] - the 
effective disorder strength in our case turns out to be larger |18| than that in the case of 
random dipoles p^,|T7|j. 

From the equations (||) it is easy to see that for g > 1/327T the fugacity y monotonously 
increases during the scaling procedure dy/dK > for all values of K. Therefore the density 
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of free charges in the system remains finite at any temperature and the KTB phase transition 
for charges never occurs. In the case of a correlated disorder @ the condition g c > 1/327T is 
equivalent to T > E c /64:ii 3 ~ 5x lO -4 ^- For typical experimental values of E c fi]J§-§[ this 
inequality can be violated only at temperatures below 1 mK. For the uncorrelated disorder 
(0) the condition g u > 1/327T includes most of physically important disorder configurations. 
E.g. if we assume that offset charges on islands are uniformly distributed between —e/2 
and e/2 we get the value g u = n 2 /3 ^> 1/327T and thus the KTB-ordered phase for charges 
should be completely destroyed by disorder. Note that this conclusion is consistent with 
the numerical results |J which demonstrate smearing of a universal jump for the nonlinear 
conductance of a 2D array with uncorrelated charge disorder. 

For g < 1/327T the RG flow (H) becomes nonmonotoneous for not very large K- < 



K < K + (see also fig.l). It was concluded [|T^] that in this case the phase diagram 
of a disordered 2D XY model shows a reentrant behavior. Recently this conclusion was 
criticized by Nattermann et al. [17] who found no reentrant behavior and, making use of the 



convention about nonrenormalizability of T (which was not used in |]T6| and here), argued 



in favour of the existence of a KTB phase under the condition equivalent to g < 1/327T. We 
believe that for our physical situation charge disorder should destroy the KTB phase even 
for g < 1/327T. Indeed in our model (in contrast to ||TB|,|T7|) the initial (unrenormalized) 
fugacity y — y is not small (yo ~ 1) and does not depend on K. Integrating the equations 
(|) it is easy to show that for g > (l/87r) exp(— 2 — 27r 2 ?/ 2 ) the line y = y can intersect 
only the RG flow lines with y renormalized to larger values (fat solid lines on fig. 1). Thus 
even if - improving the accuracy of the RG analysis - we choose somewhat smaller values 
of yo (this might only favour the KTB phase), we still can conclude that no charge-KTB 
phase transition occurs down to exponentially small g. Therefore there is a very little (if 
any) chance to provide experimental conditions for detecting the KTB phase transition for 
charges. This in fact was demonstrated in several experiments 

In order to evaluate the conductance of a 2D granular array the term St[<f>) ([I]) describing 
electron tunneling between islands should be taken into account. In the case of normal 
granular arrays this term has the form (see e.g. |T3|j2|) 

where at = n/2e 2 Rt, Rt is the junction tunneling resistance. In the weak tunneling limit 
a t <C 1 one can proceed perturbatively and expand © in powers of a t . The corresponding 
analysis is a straightforward generalization of that for a single junction Keeping only 



the first order terms and performing Gaussian integration over the phase variables (fi x for a 
given distribution of the offset charges Q x in the array we obtain 
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where F [Q(k)] = J2k Q 2 (k) / 2C A(k) is the free energy of a random distribution of charges 
Q(k) = J2x Qx ex.p(ikx) in the array for a t = 0. The conductance of a single junction in the 
array (shunted by all other junctions) can be calculated by means of a standard technique 
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l"3fl which allows to evaluate an imaginary part of the free energy F by means of an analytic 
continuation. Proceeding in much the same way as it has been done in for a single 
junction we get 



G[Q(k)} = (2R t )~ 1 exp 1 - £ 
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The expression (|TT| ) defines the conductance of a normal tunnel junction in the array with 
a fixed distribution of the offset charges Q(k). According to this result in the absence 
of the offset charges Q(k) = the conductance of each junction in the array reads G = 
(2i? t ) _1 exp(— Ec/2T). As the total conductance G to t of an ordered 2D square array of 
junctions is equal to G to t = 2G we conclude that the effective Coulomb gap of such an array 
is equal to E c /2 [[0| . 



In order to find the average value G for nonzero Q(k) one should integrate ( |TTD over all 
disorder configurations. In the case of a correlated disorder (|J) this integration is performed 
with the weight factor exp(— Fo[Q(k)]/T). As a result we obtain 



G(T) = exp[-(l/2 - lMEc/T}. (12) 

Provided the disorder in the distribution of local conductances is not very large the ar- 
ray conductance Gtot can be evaluated within the mean field approximation which yields 
G tot {T) ~ G{T) and thus 

G tot (T) « R^expi-aEc/T), (13) 

where for a correlated disorder (^) we have a = 1/2 — 1/tt ~ 0.2. This value of a is close 
to that found experimentally fl],§-§| especially if we take into account the finite size effect 
which would make the Coulomb gap somewhat larger than that for an infinite array. 

Thus we came to the conclusion that disorder in the distribution of the offset charges 
decreases the Coulomb gap of 2D tunnel junction arrays. This result has a very transparent 
physical origin. Indeed in the presence of randomly distributed offset charges Q x the effective 
charge of a half of the tunnel junctions in the array is positive. Due to that the Coulomb 
gap for electron tunneling across such junctions in one direction becomes smaller than that 
without an external charge Ec/2. At sufficiently low T these junctions yield a dominant 
contribution to the array conductance G to t for the corresponding direction of the current. 
For the current of the opposite sign another half of the junctions with the negative effective 
charge would favour electron tunneling and the Coulomb gap is again smaller than Ec- 
Hence, in the presence of a charge disorder one always has a < 1/2. The particular value of 
a may depend on the type of disorder and also on temperature (the latter is the case e.g. 
for g = g u (|)). 

The same type of arguments apply to superconducting arrays. At sufficiently low tem- 
peratures the number of quasiparticles above the superconducting gap A is exponentially 
small and we get G to t oc exp[— (aEc + A )/T}. This result is also in a good agreement with 
the experimental findings |||8|] . 

In conclusion, we demonstrated that in the presence of quenched randomly distributed 
offset charges the charge-KTB transition in 2D granular arrays and films is completely 
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destroyed already in the small disorder limit. The array conductance shows an activation- 
type temperature dependence with the effective Coulomb gap being smaller than that in the 
absence of offset charges. We believe that our results essentially explain the experimental 
data obtained for low conductance 2D junction arrays. 

We are grateful to R. Fazio and G. Schon for useful discussions and to S.E. Korshunov for 
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Fig. 1. The Rubinstein-Schraiman-Nelson RG flow in the limit of small 
disorder (g < 1/327T for the equations (8)). For sufficiently large values 
of initial fugacity y and except for exponentially small g the trajectories 
corresponding to a renormalized to zero fugacity y (thin solid lines between 
K_ and K + , K± = 7r/4 ± (ir 2 /16 — 2ir 3 g) 1 ^' 2 ) are irrelevant and y always 
increases during the scaling procedure (fat solid lines). 



